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Abstract 

We study the proximity effect between conventional superconductor 
and magnetic normal metal with a spin-orbit interaction of the Rashba 
type. Using the phenomenological Ginzburg-Landau theory and the qua- 
siclassical Eilenberger approach it is demonstrated that the Josephson 
junction with such a metal as a weak link has a special non-sinusoidal 
current-phase relation. The ground state of this junction is caracterized 
by the finite phase difference (ySp, which is proportional to the strength of 
the spin-orbit interaction and the exchange field in the normal metal. The 
proposed mechanism of the (^q junction formation gives a direct coupling 
between the superconducting current and the magnetic moment in the 
weak link. Therefore the </?q junctions open interesting perspectives for 
the superconducting spintronics. 



Usually the current-phase relations in Josephson junctions near the criti- 
cal temperature are sinusoidal j{ip) — jcSm{ip) but with lowering temperature 
the contribution of the higher harmonic terms ~ j„ sin(n(/?) can be observed. 
However if the time reversal symmetry is preserved the current-phase relation 
is always antisymmetric j{—ip) — —Jif) UJ- Without this restriction a more 
general j{ip) = josin^Lp + lpq) dependence is also possible and the generic ex- 
pression for the current in the pioneering work of Josephson |2j incorporates 
this possibility. In fact such current-phase relations have been predicted for 
Josephson coupling involving the unconventional superconductors [3] , [1] , [5] , [5] . 
The experimental verification of these predictions are still lacking. 

In the present work we demonstrate that the Josephson superconductor/normal 
metal/superconductor junctions (S/N/S) provides the realization of such un- 
usual current-phase relations j{ip) = jo sin((p -I- lpq) for the case of conventional 
superconductors when the normal layer is a non-centrosymmetric, i. e. with 
broken inversion symmetry (BIS) magnetic metal. Further on we will call this 
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junction " (pg junction" . The phase shift (Pq is proportional to the magnetic 
moment, and therefore these ipQ junctions serve example of systems with direct 
coupling between magnetic moment (internal exchange field) and superconduct- 
ing current. This opens an interesting field of application of (/?q junctions in su- 
perconducting spintronics. Varying the N layer thickness we may easily control 
the phase shift (pQ. Note that the considered situation is different from the case of 
the Josephson junction with dominating second sinusoidal harmonic - see [6], [7] 
and references cited therein. In such a case at the ground state an arbitrary 
phase drop across the junction may exist if the sign of the second harmonic is 
negative. However in these systems it is impossible to have the direct coupling 
between magnetic exchange field and superconducting phase and the properties 
of these junction are very different from that of ifQ junctions considered here. 

Before addressing the problem of the proximity effect between conventional 
superconductor and BIS magnetic metal it is useful to recall that recently BIS 
superconductors attracted a lot of attention. Namely the heavy fermion su- 
perconductor CePtaSi provides a famous example of the superconductivity and 
antiferromagnetism coexistence in the noncentrosymmetric compound [8] . Now 
the number of superconductors without inversion symmetry approaches one 
dozen and during the last years their properties were under intense studies from 
both theoretical and experimental points of view, see [5] . [TU] . pUj . . . [H] 
and references cited therein. In the presence of the magnetic field the lack of 
inversion symmetry leads to the spatially modulated helical superconducting 
phase [I0|,[I3],[I5],[I6]. 

The Josephson junctions between conventional superconductor and BIS su- 
perconductors should reveal some special features We stress that the 
aim of the present work is to study the very different situation : the Josephson 
junction between conventional superconductors with a weak link formed by BIS 
magnet. As an example of the suitable candidates for such interlayer we may 
cite MnSi and FeGe. The anomalous properties of studied junctions are related 
to the particularities of the superconducting proximity effect in the BIS metal. 

On the microscopical level the special character of the electron spectrum in 
BIS metal may be described by the Rashba-type spin-orbit coupling [18j: a{a x 
p) ■ n , where n is the unit vector along the asymmetric potential gradient and 
parameter a describes its strength. To illustrate the unusual properties of the 
BIS Josephson junction we start with a simple Ginzburg-Landau (GL) approach. 
Describing the weak link by the GL theory we assume the temperature is above 
the critical temperature of the material of the weak link and the superconducting 
order parameter is induced only by the superconducting banks. As it has been 
noted in [10], [13], the Rashba-type interaction in the presence of the field h 
acting on the electron spin leads to the following GL free energy density 
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Figure 1: Geometry of Josephson junction with BIS metal as a weak link. The 
exchange field is directed along the z-axis and the n vector is along y-axis. The 
total length of the weak link is 2L. 



F = a|V|V7 



£ n 



ft^x (^(dV')*+^*(d^))] , (1) 



where ip is the superconducting order parameter, Di = —idi — 2eAi and the 
coefficient a becomes zero at some temperature Tco : a {T — Tco). The 
special character of the BIS superconductivity is described by the last term in 
([T]) with the coefficient e ^ a. In principle the field h may be created by 
the applying external field but we suppose that BIS metal is a ferromagnet 
and h is an internal exchange field (which is assumed to be small to avoid 
the necessity to add higher order derivative terms in ([1]) [TH]) . Note that the 
origin of the spin-orbit contribution in ([T]) may be intrinsic resulting from the 
crystal symmetry or extrinsic. The latter case corresponds for example to the 
ferromagnetic layer with the in-plane magnetization and a varying thickness. 
To be more specific below we consider the intrinsic BIS metal. Schematically 
the Josephson junction is presented in Fig. 1. Further on to concentrate on the 
special properties of this junction we neglect the orbital effect. In Fig. 1 the 
magnetization is along z— axis and the demagnetization factor N=:l and then 
the internal magnetic field in the junction Hi — —AttM. Therefore the magnetic 
induction B — Hi + AttM = and then for this geometry the orbital effect is 
vanishing. Alternatively we may assume the magnetization lying in plane x — y 
(in such a case the r? vector must be along z— axis). 

In the considered case of the weak link of the length 2L , see Fig. 1, the 
order parameter depends only on the coordinate x and the corresponding GL 
equation is 

a^-,^ + 2.e/.^=0, (2) 

where we have neglected the non-linear term assuming the superconducting 
banks being at temperature slightly below their critical one. The solution of ([2]) 
is straightforward 
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■ip — A exp{qix) + B exp{q2x) , (3) 



with qi^2 = i Y ^ — e^i where £ = ^ and the condition ^ > assures 
that the weak hnk is in the normal state (i. e. the temperature is above the 
intrinsic critical temperature of BIS metal a > ac = 7£^)- To illustrate the 
particularity of the proximity effect in BIS system let us consider the contact 
of superconductor with a metal occupying the a; > half-space. In such a 
case the order parameter distribution is described by the decaying exponent 

in ([3]) V ~ exp{{ex) exp(— xy^ """" ). The difference with the usual proximity 
effect is that the order parameter decay is accompanied by the superconducting 
phase rotation. Therefore, in the weak link the phase difference proportional to 
its length would be accumulated. Due to the 2n periodicity the actual phase 
difference is limited by the interval (0, 2tt). 

To calculate the current we need to determine the coefficients A and B in 
^ from the boundary conditions at the contact with the superconductors. For 
illustration we assume that there is no barrier at the interface and we may use 
the rigid boundary conditions [IJ (i. e. the normal conductivity of the BIS 
metal is much smaller than that of the superconducting bank). Therefore the 
coefficients A and B are obtained from the continuity conditions for tp at x = 
±L: ^{±L) = |A| exp(=Fz^), with |A| being the modulus of the order parameter 
in the banks and ip is the superconducting phase difference across the junction. 
Taking into account the new expression for the superconducting current coming 
from ID) (note that there is additional contribution to the current from the spin- 
orbit term) and performing the corresponding calculation we readily find in the 

limit of the long junction ^ — >> 1 



j = 4e7|A|^y'^ - e%in((p + 2eL) exp(-2y'^ -e^L). (4) 
The current-phase relation 

j{ip)^jcsm{ip + ipa) (5) 

implies that the junction energy 

Ej jccos(¥3 + ^o), (6) 

and the minimum energy corresponds to the non-zero phase difference ip = —^q. 
Naturally the Josephson junction energy may be also directly obtained from 
the functional H]). Note that ^ describes the transition from to tt junc- 
tion when ifQ vary from to — tt. However in contrast to — tt transition in 
superconductor/fcrromagnet/superconductor (S/F/S) junctions [12], the criti- 
cal current does not vanish at the transition but remains constant. The pres- 
ence of the ground state phase difference ip^^ is a consequence of BIS and the 
h — hz component of the spin field in the weak link: </Jq — ^sM;, Note that 
we assumed the continuity of the order parameter at the bank. In general 
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the interface barrier provokes a jump of the order parameter and in our case 
it produces some additional phase rotation which may even exceed (^q for the 
large values of the interface barrier. 

The approach on the basis of GL functional thought very insightful can not 
adequately describes the systems with strong internal exchange field h > Tc, 
which is usually the case in the magnetic metals. Therefore we present also the 
theory of the (p^ junction on the basis of the quasiclassical Eilenberger equations 

m- 

Provided the spin-orbit interaction is smaller than the characteristic elec- 
tron energy scale Ep we may consider the Rashba term as and external po- 
tential in the standard scheme of derivation of the Eilenberger equations (this 
imply a << v, where w is a Fermi velocity). This approach has been suc- 
cessfully applied for the description of the superconducting state in CePtaSi 
|21j . In principle to have a complete description of the BIS junction we need 
to solve the Eilenberger equations also in the superconducting banks and take 
into account the suppression of the superconducting order parameter near the 
interfaces. The full treatment of this problem require the extended numerical 
calculations [1] . Below we would like to concentrate on the peculiar properties of 
the BIS junctions. That is why we provides the results for some cases which can 
be treated analytically. The resulting coupled equations for anomalous Green 
function fij (v, r) (matrix in spin space) are rather cumbersome for 3D case 
but strongly simplified in 2D or ID case -they are decoupled. Namely for the 
geometry in Fig.l and supposing the weak link to be quasi 2D [m x — y plane) 
the Eilenberger equations in clean limit in the region —L < x < L read 

For the junction Fig.l the superconducting order parameter in the banks may 
be considered constant as a transverse dimension of the BIS metal is small. Near 
Tc the Eilenberger equation in the bank reads (^lu + -5^^^ /12 = ^r,l, where 

is a Fermi velocity in superconductor and the equation for /21 is obtained 
by the substitute A — A. From this equation it follows that for a; > and 

> the function /12 is constant in the left bank /12 = ([A| /ui) exp(i^), while 
for < it is constant in the right bank /12 — (|A| /u) exp(i^) [22] . [23j . 

Using the corresponding continuity conditions at a; = ±L for the functions 
fij at the boundary with superconductors we may readily calculate them |22] . 
Note that the triplet components of /y vanish /n — /22 — 0. Knowing the 
Green functions permits readily to calculate the supercurrent density flowing 
trough the junction. At temperature close to the critical temperature Tc of the 
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banks at the lowest |A| approximation 

J = -ieN{0)7rT,Y, [/i2(v, (v, 

+/2i(v,x)/2+(v,a;)]), (8) 

where A''(0) is the density of state at the Fermi level. In the limit of the long 
junction L > v/h the main contribution in ^ comes from the directions \vx \ ^ w 
and the formula for the current takes a very simple form 

/ 4\h\L , 7r\ 

j M=Josm (^+^^ ^ ^ . 9 

Here jo — eiV(O)^^^ (f )^^^ ^ti*' absence of the spin-orbit interaction 

a = this expression coincides with that for j{ip) for the 2D S/F/S junction 
[21] . Comparing ^ with from GL theory ^ we see that the phase shift 
Vo = ^^T^ in both cases is proportional to strength of the spin-orbit interaction 
and the product hL. On the other hand the critical current in ^ oscillates with 
L changing its sign. This is a typical behavior inherent to the S/F/S junctions 
with the strong exchange field h >> T^. 19J. Such oscillations are absent in our 
GL approach ^ as it is adequate for h <Tc, otherwise the gradient terms in [Tj 
changes its sign and it is needed to introduce the higher derivatives terms. Such 
modified GL functional indeed qualitatively describes the oscillatory behavior 
of the superconducting order parameter at S /F proximity effect T9^ . 

In the 1 D model of the weak link (single channel approximation) the very 
similar to ^ current-phase dependence is obtained 

/ N ^w^nTtwA^ . , 4ahL, /4:\h\L\ 
.7M = e7V(0)^sin(^-|-— ^)cos(^^^j . (10) 

We considered a weak link in the framework of Eilenberger equations in 
the clean limit (ballistic regime) . In the diffusive regime the very convenient 
approach is provided by the Usadel equations |25j for the Green functions in- 
tegrated over Fermi surface -Fij(r) = (/ij(v,r)). The calculation of the current 
on the basis of Usadel approach gives us also the expression which may be pre- 
sented in the form Therefore the formation of the <pg junction by the BIS 
magnets is very general phenomenon which may be observed in both clean or 
dirty limits. 

To summarize in all approaches we obtain the Josephson junction with un- 
usual current-phase relations j{(p) = jc siTi{(p + (Pq), where the phase shift ip^ is 
determined by z-component of the internal magnetic field. Thought our model 
is applied for the weak spin-orbit interaction a << w we may expect that qual- 
itatively the phase-shift effect would be the same for the systems with strong 
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spin-orbit interaction a ^ v. In this case the characteristic length of the phase- 
shift is the same as for the — tt transition in S/F/S junctions, i. e. several 
nm. Therefore we may believe that the formation of Lp^ Josephson junction is 
inherent to all weak links/barriers with magnetic BIS metals. 

The S/F/S Josephson junction may have zero or vr phase difference in the 
ground state depending on the length of the weak link. In contrast in the ip^ 
junction the ground state is always different from zero and tt states (except the 
occasional events — mi). In the superconducting ring with tt junction the 
spontaneous current appears [26] if the parameter k — < 1, here L is the 

inductance of the system. For the ip^ junction the system energy: 

i?M = |(-cos(^ + ^o) + ^), (11) 

Therefore the minimum energy is achieved for the phase difference satisfying 
the equation 

sm{ip + (^o) + kip = 0, (12) 

which always have a non zero solution and then the (pQ junction will always 
generate the spontaneous current with the flux $ = — (f>o(f^)(l^fc)infc<<l 
limit. The SQUID with one normal and another Pq junction would reveal the 
shift of the diffraction pattern by tpQ. Note also that the (pQ Josephson junctions 
may serves as a natural phase shifter in the superconducting electronics circuits. 

The very important property of the discussed Pq junction is that it provides a 
direct mechanism of the coupling between supercurrent and magnetic moment 
- indeed the phase shift Pq is proportional to z component of the spin field. 
This means that the precessing magnetization will be directly coupled with the 
current which opens new interesting perspectives to study the coupled magnetic 
and current dynamics in Josephson junctions. Applying the voltage to the (^g 
junction we obtain the Josephson generation and the magnetic moment of the 
weak link will experience the effective field varying with Josephson frequency. 
If this frequency is close to the ferromagnetic resonance frequency, it may be 
an efficient way to generate the spin precessing. Inversely the spin precessing 
in the weak link would generate superconducting current in the circuit with p^ 
junction. 

Finally we note that even in the centosymmetric compounds the inversion 
symmetry is broken near the surface. This means that locally the Rashba 
type interaction will be present their and then the Josephson junction made by 
two superconducting electrodes attached to the surface of ferromagnetic metal 
would be a ipQ junction. 

The author is grateful to M. Houzct, E. Goldobin, J. Cayssol, F. KonschcUe 
for useful discussions and comments. This work was supported by French ANR 
project "ELEC-EPR". 
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